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RESIDUAL  THERMAL  STRESSES  IN 
GRAPHITE /PEEK  (APC-2)  LAMINATES 


By 


Y.  Jack  Weitsman  and  K.  Lee 
The  University  of  Tennessee  at  Knoxville 


Abstract 


Thermoplastic  resin  composites  undergo  a  substantial  temperature  drop  during 
their  post-manufacturing  cool-down  and  sustain  substantial  residual  stresses  due 
to  mutual  geometric  constraints  among  the  multi-directional  plies.  In  view  of 
the  time-dependent  thermomechanical  response  of  the  resin  the  residual  stresses 
exhibit  strong  dependence  on  cool-down  history.  This  paper  demonstrates  that  it 
is  possible  to  obtain  an  optimal  cool-down  path  which  minimizes  the  residual 
thermal  stresses  in  APC-2  composites. 
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1.  INTRODUCTION 


APC-2  (graphite/ PEEK)  composites  exhibit  substantial  time-dependent  stress- 
strain  response,  especially  at  high  temperatures,  and  a  very  large  disparity  be¬ 
tween  the  longitudinal  and  transverse  coefficients  of  thermal  expansion.  In  ad¬ 
dition,  the  processing  of  APC-2  involves  cool-down  from  Tj^  =  400°C  to  Tr  =  20°C 
-  which  is  about  twice  the  temperature  excursion  sustained  by  graphite/  epoxy 
composites.  The  above  factors  produce  significant  residual  thermal  stresses  in 
multi-directionally  reinforced  laminates,  which  exceed  the  range  of  linear  behav¬ 
ior.  In  view  of  the  time-dependent  response  of  PEEK,  these  stresses  are  sensitive 
to  temperature  history. 

It  seems  that  the  most  complete  characterization  to  date  of  the  time-dependent 
response  of  PEEK  and  unidirectionally  reinforced  APC-2  coupons  was  obtained 
by  Xiao  [1],  and  typical  creep  data  at  various  levels  of  temperature  are  shown  in 
Figure  I.  These,  and  supplementary  experimental  results  not  shown  here,  en¬ 
abled  Xiao  to  represent  the  time  dependent  response  of  PEEK  and  APC-2  by 
means  of  the  non-linear  viscoelastic  model  of  Schapery  [2].  However,  it  is  im¬ 
portant  to  note  that  the  creep  data  in  ref.  [1]  were  collected  under  isothermal  con¬ 
ditions  and  are  limited  to  the  temperature  range  of  20°C<T<200°C.  It  was  shown 
in  ref.  [3]  that  the  implementation  of  Schapery’s  model  requires  additional  creep 
data  under  transient  temperature  conditions.  Consequently,  Xiao’s  characteriza¬ 
tion  contains  an  uncertain  component.  The  time-dependent  response  of  APC-2 
laminates  can  be  predicted  from  the  behavior  of  a  single  ply  by  means  of  classical 
laminate  theory  [4]  upon  the  utilization  of  Schapery’s  model  and  employment  of 
the  quasi-elastic  viscoelastic  approximation  [5]. 

When  linearly  viscoelastic  materials  undergo  a  prescribed,  geometrically  con¬ 
strained  temperature  drop  over  a  finite,  predetermined,  time  interval  it  is  possi¬ 
ble  to  find  an  optimal  cool-down  path  T^ft)  which  results  in  minimal  residual 
stresses.  Such  paths  were  found  for  thermorheologically  simple  [6]  -  [8]  and 
thermorheologically  complex  [9]  viscoelastic  responses.  However,  the  foregoing 
analyses  do  not  apply  to  APC-2  composites  because  of  the  substantial  non-linear¬ 
ity  which  occurs  during  their  cool-down.  Consequently,  the  extension  of  the  op¬ 
timization  scheme  to  account  for  non-linear  behavior  is  the  subject  of  the 
present  work. 


2.  ANALYSIS 
2.1  Basic  Equations 

The  linear  thermo-elastic  response  of  a  uni-directionally  reinforced  ply  undergo¬ 
ing  a  temperature  excursion  AT  =  T  -  T]  from  some  stress-free  reference  tempera¬ 
ture  T|  (say),  is  given  by  [4] 
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El  -  ttL  Aj  =  Ol/El  -  VjL  Ot/Ex  ^  Ex  -  CxAT  =  -Vlx  Ol  /El  +  Cfx/Ex  (1 ) 

where  standard  symbols  were  used  and  subscripts  L  and  T  denote  longitudinal 
and  transverse  directions,  respjectively.  In  the  particular  case  of  a  symmetric  bal¬ 
anced  cross-ply  lay-up  a  straight  forward  employment  of  laminate  theory  yields 
the  following  expression  for  the  laminate-level  residual  thermal  stresses 

Ol  =  -Op  =  -  o: 

o  =  -  r  Ex  cxAT  (2) 

where 

1  =  0 +vxl)[i and  a=ax-aL  (3) 

Note  that  the  form  of  equation  (2)  remains  valid  for  other  symmetric  lay-ups  (in 
particular,  quasi-isotropic  lay-ups),  but  with  different  expressions  for  r.  In  the 
case  of  APC-2,  the  most  pronounced  time  dependence  occurs  in  the  transverse 
modulus  Ej.  However,  since  Ej/El  «  1,  we  may  ignore  the  time  dependence  of 
r  and  let  r  =  0.9125  to  within  ±  2%.  In  the  isothermal  case  with  the  stress  free 
temperature  denoted  by  T],  where  T  =  constant  and  AT  =  7]  -  T,  employment  of 
Schapery’s  non-linear  viscoelastic  model  modifies  equation  (2)  to  read  as  follows 
[2]: 

o(t)  =  -raV2(T)|  E(e(t)- 4(t))  [d(Vi(T)AT) /dt'd  r  (4) 

Jo 

where  the  product  V(T)  =  Vi(T)V2(T)  is  the  vertical  shift  required  to  coalesce  iso¬ 
thermal  data  (say  creep)  to  form  a  "master  curve."*  In  addition,  the  "reduced 
times",  4(t)  arid  ^(t)  in  equation  (4)  are  given  by 

^(u)=  I  dp/a(T;a(p))  (5) 

Jo 

In  equation  (5)  a(T;o)  is  the  horizontal  shift  factor  required  to  coalesce  isothermal 
data  (say  creep)  to  form  a  "master  curve."  Though  both  "V"  and  "a"  may,  in 
principle,  depend  on  o  -  it  appears  that  the  non-linearity  of  the  response  can  be 
expressed  through  a(T;o)  alone.  Under  fluctuating  temperatures,  with  T  =  T(t),  it 
is  necessary  to  re-write  equation  (4)  with  V2(T(t))  and  Vi(T(t)).  In  this  case  we 
also  have  AT  =  AT(t)  in  equation  (4)  and  a  =  a(T(p);a(p))  in  equation  (5). 


*  It  is  obvious  from  equation  (4)  that  for  T  =  constant  <j(t)  depends  only  on  V(T)  and  not  on  (T)  and 
V2(T)  separately.  This  no  longer  holds  for  fluctuating  tcm|X“ratures. 
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2.2  Reduction  of  APC-2  Response  Data 


The  isothermal  creep  data  shown  in  Figure  1  were  reduced  by  Xiao  [1]  to  expres¬ 
sions  which  fit  Schapery's  model.  However,  the  present  analysis  requires  values 
which  involve  the  relaxation  modulus  E(t),  rather  than  creep  compliance  D(t). 
E(t)  was  obtained  from  creep-response  functions  by  the  well  knowm  relation  [5]: 
1  sin(7tp(t))  ^  dloglXt) 


Values  of  E(T)  vs.  log  t  are  plotted  in  Figure  2.  The  vertical  shift  factors  V|(T)  and 

V2(T)  for  relaxation  are  taken  to  be  Vi  =  h2  ^2  =  hi’  where  h|,  h2  are  vertical 
shift  factors  for  creep.  Since  ref.  [1]  provides  only  the  product  h(T)  =  hi(T)h2(T), 

w'e  take  V(T)  =  h(T)'’.  The  horizontal  shift  factor  a(T;a)  was  approximated  by  the 
follow’ing  expression:  a(T;o)  =  ai{T}aa{o)  where  a<j(o)  =  exp  [-a  (r-'t)]  for  i  >  x  and 

a<3(o)  =  1  for  X  <  X.  a-ifT)  was  determined  from  a  set  of  values  0]=  3.3, 
log  aj,=  6.0;  62  =  2.69,  log  aj^  =  2.5;  63  =  2.1,  log  ajj  =  -8  and  64  =  1.6  log  ay^  =  -10  w'  i  t  h 

6j  =  1000/(273.1  +  Tj  °C).  Intermediate  values  were  obtained  by  linear  interpola¬ 
tion.  The  parameter  a  =  a(T)  was  obtained  by  a  spline  function  fit  of  a  (24.41)  = 
0.118,  a  (44.26)  =  0.118,  a  (74.41)  =  0.118,  a  (91.32)  =  0.146,  a  (103.1)  =  0.171  and  a 

(120)  =  0.239.  ^{T)  was  obtained  by  a  spline  function  fit  of  x  (20)  =  31.42,  x  (26)  = 
31.3,  X  (38)  =  31.03,  x  (50)  =  30.8,  x  (86.2)  =  30,  x  (108)  =  22.16,  x  (114)  =  20,  and  x  (120)  = 
17.86.  The  vertical  shift  V(T)  corresponded  to  log  V(40)  =  0.025,  log  V(129.4)  = 
-0.0031,  and  log  V(200)  =  -0.21  with  linear  interpolation  for  intermediate  values. 

Finally,  the  coefficients  of  thermal  expansion  were  taken  to  be  Ol  = 
0.5  X  loVl°C,  ttT  =  30  X  lO'^/rC  for  T  <  125°C  and  ol  =  loM°C,  aj  =  75  x  lO'^/rC, 
for  125°C  <  T  <  300°C.  To  avoid  computational  difficulties  at  T  =  125°C  w'e  as¬ 
sumed  smooth  transitions  in  values  of  at  and  aj  to  occur  over  124°C<T<126°C. 

2.3  Optimal  Cool-Dowm  in  the  Linear  Case 

The  optimization  of  cool-down  can  be  stated  as  follows:  Given  an  initial  tem¬ 
perature  T],  where  the  composite  is  stress-free,  a  final  temperature  Tp  (Tp  <  Tj),  a 
finite  time  interval  tf  and  a  thermoviscoelastic  stress-strain  relation  (e.g.  equa¬ 
tion  (4)),  find  the  "best"  time-temperature  history  T  =  Tn(t),  which  minimizes 
the  residual  thermal  stresses  o(tf).  In  all  previous  studies  [6]-[9],  which  concerned 
linear  behavior,  namely  a  =  a(T)  in  equation  (5),  it  was  found  that  Tn(t)  under¬ 
goes  sharp  drops  at  t  =  o,  from  Tj  to  T^,  and  at  tf,  from  Tn(tf)  to  Tp  , 
and  it  varies  monotonically  and  smoothly  over  the  range  o  <  t  <  tf.  Specifically 
for  the  linear  case  expressed  in  equation  (4)  (namely  a  =  a(T)  ?ta(T;o))  the  initial 
drop  is  given  by 
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(6) 


To-T,= 


Vi(To)a{To) 


a'(To)V,(To)-Vi'(To)a(To) 


which  is  a  transcendental  equation  in  the  unknown  Tq  .  Beyond  the  initial  drop, 
the  optimal  path  Tn(t)  is  governed  by  the  integro-differential  equation 


■[vi'(T-T,)+ Vi. 


dTn(t)_  E(Uf) _ 

EM  aja'[vi'(T-Ti)  + V]]-a'[vi'{T-Ti)  + Vi]  j 


(7) 


I  and  primes  denote  derivatives  with  respect 

to  the  argument.  Note  that  equation  (7)  must  be  solved  "backwards"  from  t  =  tf 
towards  t  =  o  since  E(t,tf)  cannot  be  evaluated  without  a-priori  knowledge  of  the 
solution,  with  the  exception  of  E(tf,tf)  =  E(o)  -  which  is  known.  The  procedure  is 

to  guess  Tjj  (t  f),  solve  equation  (7)  numerically  and  evaluate  T(o),  then  adjust  the 
guess  value  of  Tj^tf)  iteratively  until  we  obtain  T(o)  =  which  matches  the  root 
of  equation  (6). 

2.4  Non-Linear  Optimal  Cool-Down 

(a)  The  Three  Element  Model  as  a  Prototype  Case. 

Consider  first  the  linear  viscoelastic  c^se  with  E  (t)  =  Co  +  Do  exp 
a(T)  =  exp  (-T/A  +  B)  and  V]  =  V2  =  1.  Straightforward  manipulations  yield  [6]: 

To  =  T,  -  A,  Tn(t)  =  a[b  -  Indft)]  (8) 


In  equation  (7),  E(t,tf)  = 


if 


ds 


a{T(s)) 


where  d  (t)  =  t/X  +  exp  C,  and  C  =  B+1  -Tj/A.  The  residual  thermal  stresses  dur¬ 
ing  optimal  cool-down  are  given  by 

(l-v)o(t)/a  =  a|Co[1m  <>(0  +  Ti  /A-B]-t-Dol  (9a) 

or 


o(t)  =  -k  Tn(t)  +  k  T;  +  m  A 


(9b) 


where  k  =  aE(«»)/(l-v)'  m  =  aE(o)/(l-v).  Introduce  non-linearity  by  considering 

a  =  a(T,o)  =  exp(-T/A  +  B  -  p  o)  P>o  (10) 
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The  optimal  cooling-path  can  be  obtained  by  iteration,  using  the  linear  results  (8) 
and  (9)  as  initial  guesses  which  are  denoted  by  Namely,  in¬ 

sert  the  linear  result  (9)  into  equation  (10).  After  straightforward  manipulations 
we  obtain 


a  =  (T,o)  =  exp(-T/A”4B^’^)  (11) 

where  A^^^=  A<°V(l-pkA<°>)  ^  B^^)=  B<°> [l-p(kTi+m A'°VB^°)]  with  A^°^=  A  , 
B<o)_B.  The  primary  significance  of  expression  (11)  is  that  it  no  longer  con¬ 
tains  the  stress  o  explicitly,  hence  we  can  write  a{T,o)  =  a^^T),  and  a^^T)  is  of  the 
same  form  as  a(T).  Consequently,  (t)  and  will  retain  the  same  forms 

as  (8)  and  (9),  but  with  A^^^  and  B^^^  replacing  A  and  B.  At  this  stage,  we  substitute 
cj^^Xt)  into  equation  (10)  and  obtain  a  revised  a(T,o)  which,  as  can  be  readily  seen, 
is  of  the  form  a^-)(T)  =  exp  (-T^^Va^^*+B^^^)  where  A^^^  =  A^^^  and  B^~^=  B*°' 

[ l-P(kTi  +  mA*^’)  An  additional  substitution  yields  A'^^  =  A'^^  and  B^^’  =  B'“' 

indicating  convergence  after  three  iterations.  The  optimal  path  for  the  three  el¬ 
ement  model  with  viscoelastic  non-linearity  introduced  through  expression  (10) 
is 


Tn(t)  =  A  [b  - /«  bit)]  ,Tfi(o-)  =  Ti-A  (12) 

where  A=A/(l-pkA)  ,  B=B [l- p (kTi+m a)/b]  ,  and  b^t)=  t/A.+ exp  (i+B-Ti/a). 
The  analytic  solution  (12)  and  the  iteration  scheme  serve  as  a  model  and  a 
verification  check  for  the  numerical  scheme  which  follows. 

(b)  Optimal  Cool-Down  of  [0/90ls  APC-2  Composites. 

The  optimal  temperature  path  was  obtained  numerically,  by  means  of  a  special- 
purpose,  custom-made  computational  sub-routine.  The  procedure  was 
follows: 

Step  1:  Consider  linear  viscoelastic  behavior,  with  response  functions  E(t),a,a(T) 
and  V(T)  determined  from  reduced  experimental  data  and  with  prescribed  Tp  Tp 

and  tf.  Employ  equation  (6)  to  determine  To=T(o'^)and  solve  (7)  numerically  to 

obtain  Tolt)  ,  where  To(tf)  is  chosen  iteratively  until  a  solution  of  (7)  yields 

Tn(o")  =  To.  Denote  Tn(t)  by  fn\t). 

Step  2:  Compute  the  residual  thermal  stress  oft),  denoted  by  o^°*(t) ,  due  to  cool¬ 
down  along  TQ\t) . 
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Step  3:  Relate  o<°l(t)  to  T^\t)  to  form  c^°X=:F<o^T<°)) . 

Step  4:  Consider  the  non-linear  shift  factor  function  a(T;a)  which  corre-sponds  to 
the  reduced  data  for  APC-2.  Along  the  linear  optimal  path  this  function  is 

a(T<°);  which  can  be  expressed  as  a(i)(T) . 

Step  5:  Repeat  steps  1,  2,  3,  and  4  with  a(T)  replaced  by  a(^)(T)  to  obtain  with 

and  a(2)(T)  respectively . 

Continue  until  the  attainment  of  a  prescribed  convergence,  say 

TSjtW  +  td)  ,|e(t)|<  a  given  tolerance.  Note  the  following  observa¬ 
tions; 

1.  The  solution  of  equation  (7)  is,  in  general,  unstable  numerically.  The  optimal 
path  undergoes  rapid  drops  over  short  time-intervals,  yielding  vast  variations  in 
a(T).  Consequently,  it  is  necessary  to  employ  non-uniform  time  intervals  when 
integrating  equation  (7). 

2.  The  numerical  correlation  o  =  F(T)  and  its  incorporation  into  a(T;o) 
a(T;F(T))  alT),  as  indicated  in  steps  3  and  4  above,  require  smoothing  opera¬ 
tions  to  provide  reliable  values  for  the  derivatives  of  3(T)  to  be  used  in  equation 
(7). 

3.  We  have  no  mathematically  grounded  proof  that  the  iterated  optimization 
scheme  must  converge,  and  obviously  no  proof  of  uniqueness.  For  the  APC-2 
data  at  hand,  convergence  was  attained  after  about  4-5  iterations.  The  validity  of 
our  results  was  verified  at  least  partially  by  introducing  several  arbitrary,  small 
disturbances  in  the  optimal  path  and  comparing  o(tf) .  In  all  cases  the  resulting 
values  of  o(t,)  exceeded  the  optimal  value. 


3.  RESULTS  AND  CONCLUDING  REMARKS 

Results  for  optimal  cool-down  paths,  Tf^lt)  vs.  t ,  are  shown  in  figures  3  and  4  for 
a  cooling  time  tf=  100  min  with  initial,  stress-free,  temperatures  T|  =  250°C  and 
300°C  ,  respectively,  and  Tp  =  30°C.  These  figures  contain  results  for  the  three 
thermorheological  sub-cases:  the  "simple”  case  (V  =  V]  =  V2  =  1)  and  two 
"complex"  cases  (V  =  V],  V2  =  1  and  V  =  V2,  Vi  =  1) ,  as  well  as  three  analogous 
non-linear  cases  with  a  =  a(T;o).  Figures  5  and  6  exhibit  residual  stress  build-ups 
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during  optimal  cool-downs  which  correspond  to  the  six  sub-cases  shown  in 
Figures  3  and  4. 

Figures  3-6  demonstrate  the  important  role  played  by  the  various  material  func¬ 
tions  which  relate  the  thermoviscoelastic  behavior  of  APC-2.  Variations  among 
those  functions,  as  well  as  uncertainties  concerning  the  stress  free  temperature 
lead  to  predictions  of  residual  thermal  stress  which  vary  between  9.1  and  13.15 
ksi,  the  ultimate  transverse  stress  at  room  temperature  being  about  11.5  ksi.  It  is 
therefore  desirable  to  extend  the  experimental  characterization  work  of  Xiao  [1]  to 
include  a  higher  range  of  temperatures  as  well  as  transient  temperature  response 
to  distinguish  between  the  functions  V^fT)  and  V2(T). 
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Temperature  dependence  of  transverse  creep  compliance  of 
APC-2  unidirectional  laminate,  at  15.22Mpa.  After  Ref.  [1]. 


Compliance  (1/G 


12 


Time  (minutes) 


Time  (minutes) 
Fig.  4 


ik 


14 


Time  (minutes) 


I  I  ^  -  I  r  -  -  I  I 

tl  21  01  9  9  t  Z  0 

(]S^)  D  SSdl]S 


Time  (minutes) 
Figure  6 


